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Abstract 

This report deals with the radiation from arbitrary source 
distributions in plane stratified, anisotropic media. The formal solutions 
are obtained by an extension of modal procedures familiar from the analysis 
of isotropic waveguide regions. Special attention is given to the formulation 
of a radiation condition which requires the flow of energy away from the 
sources, and to its interpretation utilizing the refractive index surfaces 
descriptive of plane wave propagation in an anisotropic medium. These con¬ 
cepts are illustrated in detail for an ionized plasma under the influence of an 
external steady magnetic field. 
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L Introduction 

The propagation of plane electromagnetic waves in regions with 
anisotropic permittivity or permeability has received considerable atten¬ 
tion in the literature - the former in connection with the study of the 

optical properties of crystals and the electromagnetic behavior of the 

2 3 4 

ionosphere ’ , and the latter in connection with magnetized ferrites . 

In these applications, the sources of the field have been of minor concern, 
and emphasis has been placed on the determination of the wave spectrum 
in the region. However, some recent problems - for example, the in¬ 
vestigation of discontinuities in ferrite or plasma loaded waveguides, 

Cerenkov effects due to moving charged particles in a magneto-ionic 
5 

medium , or radio communication from satellites passing through the 
ionosphere — require a knowledge of the radiation characteristics of 
localized sources. 

This report deals with the radiation from arbitrary time harmonic 
sovirce distributions (varying like exp(ja)t) ) in a transversely unbounded 
region filled with a dielectric medium characterized by a tensor permittivity 
£ as shown in (1). The permeability fi^ is assumed to be constant, e may 
be a piecewise constant function of the longitudinal variable z, thereby ad¬ 
mitting linear stratification. By an extension of the modal procedure described 
for isotropic regions by Marcuvitz and Schwinger^, the transverse (to z) 
electromagnetic fields (x, y, z) and H^(x, y, z) are represented as a super¬ 
position of orthogonal transverse vector mode fields having the form V.{z)'e . (x, y) 
and L(z) H^(x, y), respectively; each mode is a solution of the source-free 
Maxwell field equations. Because of the infinite extent of the cross-sectional 
(x - y) domain and the non-variability of c with x and y, the transverse vector 
mode functions and h. have an exponential dependence of the form exp(j^x+j q y), 
where the modal index i--/(^,q) is continuously variable. As regards their 
vectorial characteristics, the modes are found to separate into two sets, *ordinary* 
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and * extraordinary* (see Appendix A), in terms of which the electromagnetic 
fields can be represented as in (7ai b). The modal amplitudes V and I are 
, found to satisfy the first-order differential equations (8) (transmission line 
equations), thereby permitting the use of network methods in their determina¬ 
tion. The continuity requirements of the tangential electric and magnetic 
fields across an interface between two different anisotropic regions, or the 
imposition of boundary conditions at a terminal surface of the region, give 
rise to coupling between the ordinary and extraordinary modal amplitudes. 

Their systematic calculation via equivalent network procedures is emphasized 
in Sec. n. 

The transmission line equations (8) descriptive of the spatial behavior 
of y and I in each layer contain as parameters the propagation constant 3^ and 
the characteristic impedance Z. As shown in Appendix A, 3^ in an anisotropic 
medium is generally a complicated multivalued function of the transverse wave- 
numbers (§, q). Its proper definition on the integration path, essential for a 
unique specification of the fields via the modal representations (7a, b), is ac¬ 
complished by recourse to a radiation condition which requires that the energy 

7 

flow due to a localized source distribution is outward from the source region. 
The associated restrictions on3(, complicated by the fact that the directions 
of energy and phase propagation generally differ in an anisotropic medium, are 
discussed in Sec. Ill - - both analytically, and through use of the refractive 
index surfaces for the medium which specify the variation in refractive index 
as a function of the direction of propagation of a plane wave. Special attention 
is given to the case of a plasma under the influence of a longitudinal d. c. 
magnetic field, and to the effect of the plasma parameters in determining the 
propagation characteristics of electromagnetic waves. The definition of the 
multivalued function for various plasma parameters is examined in detail 
in Sec. IV, thereby rendering unique the formal solution of this class of 
radiation problems. 

The results of this analysis have been applied elsewhere”* ~ to the 
detailed study of the radiation field of an electric current element in an infinite, 
and semi-infinite, homogeneous, anisotropic plasma medium. 

Investigations by other authors have been concerned primarily with 
the study of radiation from infinite, semi-infinite, and single-slab anisotropic 
plasma regions. For orientations of the gyrotropic axis perpendicular to the 
interface, the formal solutions for these configurations emerge as special cases 
of our analysis. While the imposition of a radiation condition is essential in the 
formulation of radiation problems in infinite and semi-infinite regions, the de¬ 
tailed discussion of its effect on the proper specification of the representati^on 
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Integrals as carried out herein has generally been omitted by other authors. With 
reference to the infinite medium problem we cite the work of Abraham^^, Bunkiri^\ 
Kogelnlk^^, Kuehl^^, Mittra^^, Mittra and Deschamps^^, Clemmow^^* Motz and 
Kogelnik^^, Chow^®. Abraham* has employed a two-dimensional Fourier integral 
representation for the fields similar to that presented here, and has also pointed 
out the utility of the refractive index surfaces in dealing with the radiation condition; 
his analysis in connection with the latter does not; however, enter into the function- 
theoretic questions treated in Secs. Ill and IV. The approach of the other authors 
(except Clemmpw) differs In that they proceed via a three-dimensional Fourier 
integral formulation which, though imbued with a certain formal elegance for the 
infinite medium problem, appears less convenient for an asymptotic analysis of 
the radiation field, and also seems not directly suited to the study of stratified 
media. Radiation from dipole sources in the presence of an anisotropic plasma 
half-space has been analyzed by Barsukov Arbel , and Tyras, Ishimaru, and 
Swarm Hodara^^ has recently been concerned with the problem of radiation 
through an anisotropic plasma slab; he employs a two-dimensional Fourier trans¬ 
form procedure but omits any function-theoretic discussion of the type mentioned 
above. It also seems that the matching of the boundary conditions at the slab inter¬ 
face and at the source could be more easily achieved by the network procedure 
described here, as may be seen from the analysis of Wu who has used the network 
approach for the slab problem. Several studies of the simpler problems of radiation 

^ O A 9 C 

by specially oriented line sources in a gyrotropic medium ’ ', or of radiation in 

a plasma subjected to weak^®{u>^“* 0 in (35)) or stron^^®* w)external 

magnetic field have likewise been carried out. 

Green's function representations for general anisotropic waveguides 
have been obtained by Bresler and Marcuvitz via abstract operator methods. For 
the special class of problems herein, our procedure, though less general, yields 
the result more directly and by conventional methods. 


The authors are indebted to G. Meltz of the Air Force Cambridge Research 
Laboratories for calling this reference to their attention during the preparation 
of this manuscript. 
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a. Reduction of the field equations 

Consider an arbitrary (but prescribed) distribution of time-harmonic 
sources of electric current J(£) and magnetic current M(r) in a medium com¬ 
prising a series of transversely unbounded^ parallel anisotropic layers (Fig. 1). 
A coordinate system is chosen so that the z-axis is perpendicular to the layer 


( See Fig. 1 in Appendix C) 


interfaces! and each layer is assumed to be characterized by a scalar 


permeability u and a tensor dielectric constant 

?<•>= .*<•> + > z , .*<“4/ 

» -t — 0—0 z -t “I . . 


2 


where ^ ^ denotes the a-th layert it ^ is a transverse (to z) dyadic whose 

representative in an x-y coordinate system is as shown, and z is a \mit 

— o 

vector along the z-direction. These dielectric tensors have the form appropriate 
to an ionized plasma medium under the influence of a d. c. magnetic field along 
the z-axis. We seek a solution of the steady-state Maxwell field equations 


V X - M<“> , V X 


in each layer, subject to the required continuity of the transverse electromagnetic 
fields E^^^and ^ at each interface, to prescribed boundary conditions at the 
z-termini of the region (if any), and to a radiation condition at infinity. A time 
variation exp(j(,>t) is understood throughput. 
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If the configuration in Fig. 1 is viewed as a waveguide with axis 
along z, the formal solution of the boundary value problem in (2) can be 
effected by an application of guided wave techniques developed for isotropic 
and anisotropic regionsFirst, by taking the scalar and vector products of 
Eqs. (2) and the longitudinal unit vector one derives after some rearrange¬ 
ment the equations for the transverse field components (note: 1* = e^ 1), 


-TT = j“^'o it + 


^«t 

dz 


^t + 


V V 

^ t^ t 


V V 

t t 


H X z + M X z , 


— t —o 


-t -o 


z xE. + z xj. 
— o —t — o — t 


from which the longitudinal components can be obtained via the relation, 


E = 
z 


;i- V-(Hxzj-jl, H= ~ V • ( z X E J -M I . 
L t \-t -o/ zj z I t ^_o -t) zj 


To simplify the notation, the superscripts identifying the various layers 
have been omitted, and the following definitions have been introduced: 


T it*^ 


e = c = e c 
z o z 


1 2 ,2 
k = k e 


1,2 2., 

^o ' “*^o®o ’ 


where is the dielectric constant for free space, c is the normalized 

longitudinal dielectric constant c /i , and c, _ = f, _/c are normalized to 

2 0 1, Z 1, 2' 

= (V - d/dz) represents the transverse gradient operator, 1^ the 
transverse unit dyadic, and 


It -It 


T— z X V. M 
jwp, — o t z 


M=M, + ^z xV, J , 

— t —t JWE—O tZ 


are equivalent transverse electric and magnetic source current distributions. 

Next, we seek a representation of the transverse field components in 
terms of a set of transverse vector eigenfunctions which are complete, orthogonal, 
and individually satisfy the homogeneous field equations and the required boundary 
conditions, in the transverse domain. As shown in Appendix A, the electric mode 

For the plasma case, normalization to c implies that c,. = F in the absence 
of an external d. c. magnetic field {see Sec. IV). 
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functions "e and the magnetic mode functions K comprise a continuous dis¬ 
tribution of plane waves whose vectorial characteristics can be grouped into 
two categories, "ordinary" and "extraordinary", to be denoted by the subscripts 
o and e, respectively: 


. ^-jk(g X + ti y) 

= ®o.e<^*’'>. -Zi- 


" , (4a) 


5o,e<*’ 


h 

— o, e 


j.--jMS X + t»y) 


(4b) 


where 


e 


+ A 
+ A yfz 


n + 

- 2 - + -"ftr 




£ = ^ 


^o * ®o ■ 7=^ ^ ^ A ' " — » (4d) 

n a ~ 2 V2 ‘ 


e 

and 

6 = 


bo 


‘2 A _ r 4 . ^,2 7 , 2 

r-r * A = v/o + 46 (1 - CT 

'l " ^ 


) 8gn(c, -.1). 


(4e) 


While the eigenvalue problem could have been treated for general dissipative 
30b 2 2 

media, the above analysis has been restricted to the lossless case for which 
c, and C 2 are real. It is noted that the wavenumbers ^ and t) have been 
normalized to k, thereby implying that e > 0 and k> 0. The modifications 
introduced when c < 0 are discussed at the end of Appendix A. The orthogonality 

"properties of the position vectors e and h differ according to whether A 
is real or imaginary (the square root in (4e) is defined to be positive when real 
and negative imaginary otherwise): 


e • h X z = 1 , 
—o —o —o 
e e 


e • h X z = 0 , 
—o —o —o ’ 

e e 


e • h X z =0 

-g "S 


h X z =1 
-8 


A real 


A imaginary , 


(5) 


where the asterisk denotes the complex conjugate. Thus, the mode functions 
in {4a, b) satisfy the bi-orthcqonalitv relation. 
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/ dx /dy e*(x, 5 (x»y;4'. n') X = e*. h x z ^ 6(^-g‘)8 {t|- tj') , (6) 

-00 -oo 

which, upon insertion of the appropriate modal subscripts o or e, is simplified 
further from (5). 

The various vector functions in (3a, b) can now be represented as 

follows: 


E^(x,y,z) = /d^/dtj |^V^(z;^,ti)£^(^,Ti) + Vg(z;4,Ti)e^(g,T|)^ 


.-jk(ex+t,y) 

(7a) 


[ '1 


■jk(ex+t,y) 

(7b) 


z X J (x, y, z) = -^ fd^Jdn H (z;^ , Tl)h , H) + i^(z;g , t^)h , t)) I 

-00 -00 L 4 (7 c) 

Mj(x, y, z)x£^ = /d4/diirv^(z;g,il)£^(^,7l) + v^(z;^ , tj) e^(^ , Tj)^e"j^ 


-jk(^x+ny) 

(7d) 


The completeness and existence of these representations (for spatially confined 
source distributions) in the (x-y) function space follows from the theory of the 
Fourier integral and from the field behavior at infinity, while that in the 2x2 
vector space is assured by the linear combination of the e- and o- eigenvectors. 
Upon substituting Eqs. (7a-d) into (3a, b), interchanging the orders of 
differentiation and integration and noting that the operator can then be 
replaced by -jka, one may equate the resulting Fourier transforms on both 
sides of these equations. Dot product multiplication of the first and second 
of these equations with (h x z ) and (z x e. ), respectively, and use of 

Oj C O ““O •*0| 0 

Eqs. ( A4 ) in Appendix A and of the orthogonality relations (5) then yields the 
following expressions to be satisfied by the transforms V and I : 

Of 0 ^ 


- -3f~= jkJfZI + v , - = jk3?YV+i 


( 8 ) 


These equations, valid separately for both the o- and e- modes (for real or 
imaginary A), can evidently be interpreted as transmission line equations, 
with V, I, Z, and Y playing the role of voltage, current^ normalized 
propagation constant, characteristic impedance and characteristic admittance. 
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respectively, while v and i represent voltage and current generator 
distributions (Fig. 1(b)). and Z = 1/Y are determined from the eigenvalue 
problem in Appendix A as 


<"o “ V 

e 


where 
U 


fUl j 


- W 


^1 ■ 


^1 2 

- a , 


Y = - 14 -- 

e e J: a A) 


w = £,(1 - - a ^) . 



S = 


^1 ■ ^2 


, {8a) 


f~z -' ^ 

V U - W can alternatively be written as - l)(A/2). The voltage and 
current generator strengths are evaluated by inverting Eqs. (7c, d) via (6): 


. (8b) 


00 00 


00 00 


«00 -00 


"" f ?o * ^t ’ ^o ^ /dx/dy h o • real; 

* -oo-oci e 


(9) 




For imaginary A, e and h are replaced by e and Ti , respectively. 

— o, e —o, e e, o — e, o ^ ^ 

From the two-dimensional divergence theorem 




00 00 00 00 

JdxJ^dy_A"V f= - J dxf dy f V ’ A + ^ ^ H • 

-00 -00 ^ -00 -00 *■ s 

where A and f are suitably continuous vector and scalar functions, s denotes a 

contour bounding the transverse cross-section at infinity in the x-y plane and n 

is a unit vector normal to s, one observes that if f = 0 on s, the double integral 

on the left-hand side is equal to the double integral on the right-hand side. Upon 

applying this result to Eq. (9) (after substituting (3e) and assuming that J and M 

2 2 

are spatially confined sources), one obtains alternatively, 


1 

o 

e 


O') 00 

^ I f It' 

-00 - 00 ^ 


— 

e 

— o 
e 


dx dy 



— Jr 

M e 
z —o 
e 


a X dx dy 


(9a) 


'^o = / / 'iig - 7:^ 7 /Jzlo^^o'-? • ( 9 b) 

e -00 - 00 6 -00 _oo “ 

which result can now be applied also to discontinuous current distributions. 

The proper definition of the multivalued function “if ( 0 ), and the dis¬ 
position of the integration path in (7a, b) in relation to its singularities, will 

One observes from Eqs. (4c, d) that the analytic continuation of the functions 
e* and h* from the domain of real A into the domain of imaginary A yields the 

functions ^ and h* respectively, evaluated for imaginary A. Thus, it suf¬ 
fices to evaltiate the alscontinuously represented inner products in (9) for real A 
and continue these functions into the imaginary-A domain. The same is true for 
other inner products encountered later on. 
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be discussed further in Sec. Ill, For the present we note only that9fo and 
are associate^ respectively, with the + and - signs in (8a), and that t/u^ - W 
is defined to be positive when real. When imaginary part, the 

associated wave is non-propagating (along z), while real represents a 
propagating wave. Since i and v are specified in terms of known quantities, 
the solution of the Maxwell field equations (2) has been reduced to the solution 
of the transmission line equations (8) in each layer, subject to the required 
continuity conditions at each interface and to the specified boundary conditions 
at the z-termini of the region. 

b. Solution of the modal network problem 


At an interface between two different anisotropic regions, the re¬ 
quired continuity of E ^ and in (7a, b) can be achieved by assuring the 
continuity of their Fourier transforms. If parameters pertaining to the two 
regions are distinguished by superscripts a and P, respectively, the boundary 
conditions are satisfied if 


kOOLW eW . vW e = k(»> [v'W . V<W .<« ] , 

,(W [,(.),(«) , ,(a) ,(a,j , ,,a) J,,,) ,,p) ^ ,(p, J _ 


(10a) 

(10b) 


provided that the normalized transverse wave numbers in the two regions are 
related via 

kWjW , . * (10c) 


In view of the orthogonality relation (5), Eqs. (10a, b) can be reduced to the 
matrix form: 




v(P) 


t t „ 

T(a) 
L ^ J 

= T 

z 

, 

^p) 

.. 


ao ap 

- ‘ep ‘pp . 


where V and~I are the column vectors 

v<‘' 

r(i) 


v<‘>- 


J 


T<‘>. 


l(i) 

O 


P , 


( 11 ) 


(11a) 


while T is the impedance transfer matrix. The matrix elements t., , 

. . ^ ij ’ 

1» j ~ a, p, are themselves 2x2 matrices whose elements are given by: 


* It is implied that q in the integral representations (7a, b) are replaced 

throughout by the appropriate i = a, p, whence = V (z: 

- V / t(P) (3) o, e o, e 

o, e " o, e'^’ § » 'H ), etc. This is to be borne in mind when interpreting 

formulas involving parameters with different superscripts. 
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*op = 


0 , 


(lib) 



r (p) 

,^(a)* (P) " 

h' ' X z • e'^' 

M O ** O 



V (a)* (p) 

h' ' x z • e'*^' 

— e —o —o 

1.(4* . (P) 

h' ' X z • e''^' 

— e —o —e _ 

, (11c) 


k<«i 

(a)* t-(P) 

e' ' • h ' X z 

(a)* .(P) 

lo • * ^o 


‘pp"^‘ 

(a)* .(P) 

e' ' • h X z 

— e —o — o 

(a)* .(P) 

e' ' • h ' X z 

—e — e — o 

(lid) 


Expressions (11c, d) apply for real or imaginary ^ (see footnote on p. 8). Since 
the submatrices t^^ and t^^ are not diagonal, the ordinary and extraordinary 
modes are coupled at the interface. Because t^^ = t^^ = 0, the coupling occurs 
in a particularly simple manner and can be schematized in terms of the trans¬ 
former network shown in Fig. 2. It is also noted, in view of the symmetrical 
character of Eqs. (10) as regards a and P, that the elements of the inverse 
matrix T"^ are given as in (11c, d) provided that the superscripts a and p are 
interchanged throughout. 

In the interior of each slab region, the equivalent modal network 
comprises the two transmission lines representative of the o- and e-modes, 
respectively, as shown in Fig. 1(b). If the slab in question has a length d, 
the impedance transfer matrix for the slab region is obtained from simple 
transmission line theory as follows: 


*^12 

‘zi ‘Z2 


( 12 ) 


where the 2x2 submatrices t.j, i, j = 1, 2, are diagonal and are given by 
tjj ~ ~ cos (k*)^ d) , Y t 22 ” Z t^^ ~ j sin (k d). 


(12 a) 


^ and Zare the diagonal propagation constant and characteristic impedance 
matrices, respectively. 



(12b) 
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The matrices in (12a) are interpreted by recalling that if^ is a diagonal 
matrix, with elements a matrix f()f) is also diagonal and has as its 

elements f (V .) . The voltage and current vectors and at the slab 

face z = are then related to the analogous quantities at z = z^ via 


“ yd) ’ 
- - 


= T 


V(2) 

T (2) 


^2 " h 


Repeated application of Eqs. (11) and (12) allows one to express the voltages 
and currents at any point z^ in the region in terms of the voltages and 
currents at any other point . The overall transfer matrix descriptive 
of the network between and z^ is then composed of the ordered product 
of the transfer matrices of the network constituents in this region. 

If the region is terminated at z = z^ in a plane surface on which 
and are related by the boundary condition 


E. (PI z ) = Z • H. (p, z ) x z , I 

— t o — — t — — o — o 

where the transverse dyadic Z denotes a constant anisotropic surface im¬ 
pedance whose representative in an x-y coordinate space is 


Z -9 i 


then the corresponding relation between ^ and ^ at z = 0 is given 


via (7a, b) by 


V e + V e = Z . (I h + I h ) X z 
0—0 e— e — o — o e — e —o 


Use of the orthogonality relations (5) allows one to deduce the terminal 

A. 

impedance matrix Z ^ for this structure as 


where for real or imaginary A (see footnote p. 8), 

r * * 

h xz *Z»hxz h xz*Z*hxz 

— o —o - —o —o —o —o — —e —o 


Z 

s 


. * 
h X 
— e 


z • Z . h X z 
-o ~ iio -o 


h X z • Z • h X z 
— e —o ~ —e —o 


* In this equation, the superscripts refer to quantities at different z - locations 
in the same region. 



f! 
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It is to be noted that even a scalar surface impedance Z (for which Z = 1. Z ,Z 
constant) couples the o-and e-modes (see Eqs. (4ctd)). Only when Z^ = 0, » 
(short or open circuit) or when (7= 0 (normally incident plane wave) does the 
coupling disappear. ^ 

It is frequently more convenient to deal with a traveling wave representa 
tion involving incident and reflected waves rather than with the standing wave re¬ 
presentation employed above. The required transformations* well-known in 

3Z 

linear network theory, are summarized below. Let 

V = a + E , "I = Y (a - S) , (16) 

where the wave vectors a and ¥ distinguish the amplitudes of waves traveling 
to the right ai^ left, respectively (see Fig, 3), 

. (16a) 


and Y is the characteristic admittance matrix (see (12b)'). Then the scattering 
transfer matrix T^ provides the direct connection between the incident and 
reflected wave amplitudes at terminals a and P as follows 


:(a) 


-(a) 




X(P) 


aa ap 

’■pa ""PP 

where the scattering transfer matrix T^ is related to the impedance transfer 


(17) 


matrix T in (11) via the linear transformation 
s 


1 

A 

1 

1 

1 

A 

1 

= “2“ 

A 

1 

m 

* L -Y<p) 

^(P) 


(18) 


1 denotes the 2x2 unit matrix. When applied to an interface between two 
anisotropic media, t^p = tp^ = 0 from (11b), and (18) yields the following 
expressions for the 2x2 submatrices 


aa 


'aP 




t + Z^®^ t-ftY^P^ 

aa pp 


aa 


- Z<°> t Y 

- Z tpp Y 




(19a) 

(19b) 
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For a slab of length d, the scattering transfer matrix T corresponding to 

A ® 

in (12) has the simple diagonal representation 




^22 


® » 


( 20 ) 


A 

with defined in (12b). Like the impedance transfer representation, the 
scattering transfer formulation is well suited to the analysis of cascaded 
networks as- encovmtered in stratified media. 

If the sources are located in a semi-infinite medium, the reflection 
phenomena are analyzed conveniently in terms of a scattering matrix representa¬ 
tion which expresses the outgoing waves at terminals (a, p) in terms of the 
incoming waves via the scattering matrix ^: 




-f 


ti-) 

^P) 


8 S a 

aa ap 

■pa "pp J 


( 21 ) 


The scattering matrix elements are related to those of the T matrix via 


-1 

■aa ' “^ap "^pp 


■ap ~ ^aa ” ^ap ^ PP ^pa 




■^Pa 


.(2 


The above equations can also be used to describe an interface between 
an anisotropic and an isotropic medium. In this instance, it is more convenient 
to employ in the isotropic medium the linearly polarized E and H modes described 
in Appendix A. 


(See Fig. 3 in Appendix C) 


■(■Explicit expressions for the various coupling matrix elements, obtained 
after carrying out the operations and substitutions in the text, are given in 
reference 8. These rather cumbersome formulas are not listed here because 
of space limitations. 
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ni. Radiation condition and specification of ^ 

It was noted in connection with Eq. (8a) that is a multivalued 

2 2 1/2 ® 

function of + 11 )' = 0 ; and that its analytic properties must be speci¬ 

fied if the integral representations (7a, b) for the fields are to be rendered 

if 

unique. The definitions 



\/u + (ej - 1) A/2 , 

(23a) 

5 yu - Vu^ - w « 

U - (ej - 1) A/2 , 

(23b) 


with U, W, and A specified in (8b) and (4e), respectively, imply that has 

2 ” 
branch point singularities at U = W (A=0) and at those values of a for which 

X * C, while TH. has branch point singularities at A = 0 and at those values 

O 6 

of O for which ^ - 0, Although the branch point at A = 0 occurs in both the 

® ifif 

ordinary and extraordinary integrals in (7a, b) , the sum of the ordinary 
and extrordinary integrands is an even function of A and therefore single¬ 
valued at A 3 0. This evenness, resulting from the fact that all ordinary 
quantities differ from the corresponding extraordinary ones oxily by the 
algebraic sign of A (see {4c, d), (8a, b) and (9)), implies that a series ex¬ 
pansion of the combined integrand about A s 0 contains only even powers of 
A, whence A s 0 is a regular point. Thus, no special care need be taken in 
the definition of ^ at A = 0 as regards the total integrand; if the ordinary 
and extraordinary integrals are treated separately, one may introduce con¬ 
venient branch cut configurations relative to the branch points at A = 0 which 
render the integrands single-valued on a Riemann surface associated with 
this sing’Uarity. 


In plane stratified media, V and I may also possess pole singularities 
on the integration path. For the case of a single interface, these (surface 
wave) poles have been discussed in reference 8, 

♦♦The two terms in the integrand of (7a) give rise to separate contributions 
which will be called ordinary (subscript o) and extraordinary (subscript e), 
respectively. Analogous considerations apply to (7b). 
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For a region comprising a series of slabs of finite width along the 

z-direction, the voltage and current solutions may be even functions of the 

propagation constants in each slab. In this case, the points 

where X - 0 and X'®' = 0, are regular. The situation is different, how- 

ever, if one of the regions extends to z = ±. «. In this instance, the integrands 

are not even functions of the appropriate propagation constant, and the matter 

of the definition of X must be studied in detail. The investigation is di- 

o, e ^ 

rectly connected with the specification of the boundary condition at | z|— 

Consider a homogeneous anisotropic region which occupies the half 


space < z < <», and assume that all sources are located in the space 

z < Z 2 where Z 2 > Then the solution of the transmission line equations (8) 

at points z > z^ comprises traveling waves characterized by the functions 

exp (-j z), where k and have the values appropriate to the region 

0|© .0^© ^ 

in question, and k is assumed to be positive. The integration in (7a, b) extends 


over all real values of a (see Appendix B). One notes from (23a, b) that 9^ 
and X may be either real, imaginary, or complex, depending on whether 

2 ® 7 

JCq g i® positive, negative, or complex (the latter case obtains when < W, 
i.e. , A is imaginary). For those values of o for which or jfg has a non¬ 
vanishing imaginary part, the requirement 


assures that the fields remain finite as z<-4 00 and serves to define the multi¬ 
valued functions ^ The associated mode fields decay exponentially with 

increasing z and represent a non-propagating wave. 

If or is real, the exp{-jk ^z) solution represents a pro¬ 
pagating plane wave. Since all sources are confined to the region z< z^, there 
must be a net flow of power toward z = «> and it appears plausible to suppose 
that the power flow vector for each constituent propagating plane wave likewise 
has a component along the +z direction. It will now be shown that the latter 
condition follows from the former, i. e., that net total power flow toward z = « 
implies that each propagating plane wave carries power in this direction. This 
radiation condition (see also reference 7) then permits the unique determination 
of K. g over the range of a for which the propagation constants are real. It 
must be emphasized that the "energy radiation condition* is distinct from a 
"phase radiation condition" since the directions of phase and energy propagation 
in an anisotropic medium are generally different. 

i-7“^-:- 

In view of the remarks at the end of Appendix A, the conclusions derived herein 

can easily be adapted to imaginary k. In that case, exp(-jk9Cz)—♦exp(-j|k|^<z), 

etc. 
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The requirement of net outward average power flow through a 
plane at z = c> z^ t where c is a real constant, can be phrased as 

°° °r * 

P = Re / dx J dy £ • S(x, y, c)> 0 , S = E x H , (25) 

^ .00 - 00 t i 


where S is the complex Poynting vector. Upon substituting the integral 
representations {7a, b) into (25), interchanging the orders of integration 
and recalling the orthogonality relation (6), one obtains: 


P = Re/ ffdi dr, [' 
real 


I + V* I 1 + r /d^ dr, f V* I + V*I 1 > > 0 , 

o e ej 'toe e oji -• ’ 

Tf, conrplex j 


(26) 


where the first integral extends over that portion of the - r, pl4ne for which 
2 

is real (A real), while the second extends over the remaining region 
o, e 2 

wherein is complex (A imaginary). Since 

O I c 


V 

O, 



c 

o. 



z>Z2 


(27a) 


where C and C are constants, one notes from the transmission line 
o e 

equations (8) (or from the equivalent network picture in Fig. 1(b) ), that 


I 

o. 



Y 

o. 


V 

e o. 



z>Z2 


{27b) 


so that (26) can be written as: 


= ■‘I IfolVol^ + YjvJ^j + //d{ df, 

[iJCreal ■j^comiiex 


Y V V + Y V 
e e o O' 


(28)^ 


From (23a, b), 

2 *2 

Xq ’ when A is imaginary, (29a) 

Hfi 

so that Xq - i Pfg • Since and must both have negative imaginary 
parts (from (24) ), it follows that 




-x: 


Im X < 
o, e 


>0 . 


when A is imaginary. 


(29b) 
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whence, from (8a), = - Y^. Hence, the integral over the range of 

complex , which arises from a coupling of the ordinary and extra¬ 
ordinary mode energies, is imaginary and does not contribute to P . If i 
and (or) is negative real, i. e., and (or) Jf. i® imaginary, 
then Y^ and (or) Y^ is also imaginary (see (8a) ) and the resulting integral 
does not contribute to P^. Hence, (28) can be simplified to: 

P = ff ^ol^ol^ de dti + ff YjVgl^ dg dr, >0 , (30) 

jC^real real 

where the equality has been omitted since we are considering propagating 
waves which carry energy along the +z - direction. The symbol "Re" is 
superfluous because the integrals are now real. Since the radiation con¬ 
dition (30) must be valid for arbitrary source distributions iii the region 
2 < z_, (i. e. , for arbitrary V and V ), each of the integrands must satisfy 
the inequality whence P^ > 0 if, and only if, 

Y >0 when is real. (31) 

0| 6 ' Of 6 • 

Hence, as stated above, each propagating plane wave must individually 
satisfy the radiation condition. 

(8a) can now be employed to determine the algebraic sign of 
^ in (23a, b) when ^ is real. Since the characteristic admittances 
must be positive, the sign of ^ is identical with that of (A^ + a^A). 

This condition evidently depends on the values of the constitutive para¬ 
meters , € 2 * and is investigated below in detail for a plasma medium. 
Eqs. (24) and (31) specify ^ uniquely for all real values of g and q 
(or a), and can be used to determine the analytic continuation of the 
function ^(a) around branch point singularities located on the real o 
axis. Some general remarks can be made concerning the determination 
of the range of propagating modes for which 

Vo = U± - W (32) 

e 

2 

is positive real. One evi dently must have U - W>0. IfW>0 for 
some value of a, then x/u^ - W < |ul and < 0 if U ^ 0; thus, 

both the ordinary and extraordinary modes propa gate whe n U > 0, and 
neither propagates when U< 0. If W < 0, then Ju^ - W > (u|; in this 
instance, the ordinary mode propagates for U ^ 0 while the extraordixiary 
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2 

mode does not. The points ^ at which U - W = 0, or A = 0, are obtained 
from (4e) as 



(33) 


If 6 > 1, these points lie on the real a axis and give rise to a- intervals 

for which is complex. 

0| o 

While the specification of the function ^ (o), when real, can 

be carried out analytically from a study of (31), there also exists a simple 

graphical procedure utilizing the dispersion curves 9((o) vs. a which locate 

the real solutions of the dispersion relation in the o plane (see (A9) and (A15) ). 

For a plasma m.edium under the influence of a steady external magnetic field, 

these dispersion curves fall into various separate categories (distinguished 

by different ranges of the values of the applied signal, plasma, and cyclotron 

frequencies) and have been investigated in great detail in connection with plane 

2 3 33 34 

wave propagation in the ionospherfe.' In particular, it has been shown that a 
plane wave characterized by the variation 


- jk 

e •’ 

where £ = x^ i ■*’Zo ^ ^ —o ^ wave normal, carries real power in 

a direction perpendicular to the dispersion curve at the point (^ , ri, ), 
and that the angle between the real part of the complex Poynting vector S 
and the wave normal is less than 90°, e direction of energy flow is 

commonly called the "ray** direction, and the above-mentioned relations 
are schematized in Fig. 4 for a typical case for which the dispersion curve 
has a closed and an open branch. The dispersion curves (or refractive 

3t< 

index curves ) can be shown to be rotationally symmetric about the d. c. 


i x + q y+lK(a)z 


= e 


- jk£. r 


(34) 


♦Upon writing exp j k£ * £ ( = exp ” j P k • £ J , where k is the wave- 
number defined in (3d) and the vector k is parallel to £, one may interpret 
p as the refractive index (with respect to a medium having a dielectric 
constant e = e e ) for the wave traveling in the direction k. Since p = X + o , 
the distance from the origin to the X vs. G surface yields the magnitude of p 
as a function of the polar angle 6 measured from the direction of the applied 
d. c. magnetic field. Hence, this surface also constitutes the "refractive 
index surface" for the medium; the latter designation is used extensively 
in ionospheric propagation theory. 
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I 


(See Fig. 4 in Appendix C) 


magnetic field direction (z axis) whence a cross-section in the x-z plane 
suffices. In Fig. 4, the spatial coordinate axes have been superposed 
upon the 3^ - ^ axes, so that the vector directions of the wave normal 
and ray are directly those in the x-z plane. 

The radiation condition (31) can now be interpreted as cor¬ 
responding to those values of “X- and a for which Re S has a component 
in the + z direction. Since £ • (Re S)> 0, the pertinent segments of the 
dispersion curves are those shown shaded in Fig. 4. One notes that it 
is possible to have £ • = X < 0 while (Re S) ' 0i corresponding 

to a "backward" wave in which the directions of energy and phase propagation 
along z are opposite. While this simple graphical construction allows the 
identification of those portions of the dispersion curve which contribute 
propagating waves carrying power in the positive z-direction, it gives no 
information as to which curve segments correspond to Xq or Jfg defined 
in (23a, b). The assignment of the proper ordinary and extraordinary 
branches defined herein requires a further study of the behavior of the 
functions U and W as discussed in connection with (32). These remarks 
are developed further in the next section where we analyze in detail the 
nature of the dispersion curves, and the associated disposition of ^ * 
for a plasma medium. 
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i 

i 

5 


rV. The function y(a) for a plasma medium 

An ionized plasma medium under the influence of a steady external 

magnetic field H along the z-axis can be characterized by the following 
* ■ ® 2 3 

constitutive parameters: ^ 

2 

u (i> u 
C ■ 


2 2 
U 0) 


1 


1 + 


C p 

77—ZTT—z: 

(w^ - « )(u - Wp) 


7- 2 —7-2~ 

(a>; - 


= 1 


—^ * ( 35 ) 


where w, « and w are the angular applied, plasma, and cyclotron frequencies, 
P c 

tespectively. In this simple description of a plasma, only the electrons are 
considered mobile, and collisions with ions and neutral particles are neglected; 
also, the impressed a. c. field amplitudes are required to be small. In terms 
of the electron density N, the electronic charge e, and the electronic maei *n, 
one may express Up and as: 


OJ = 


eN 

me. 


OJ = 


£ u H 
^o o 

m 


(35^) 


To assess the properties o{^a.» a function of o, it is necessary to 

investigate the behavior of U and W (see (8b) for t 0) for various values of 

u, «p. To facilitate this analysis, the parameters c^, 2ej/(ej + 1), c^, a^, 

have been plotted in Fig. 5 as functions of «. Fig. 5(a) pertaing to the case 

« < w , while « > w in Fig. 5(b). If a scale for o^ is superimposed along the 
P ^ P ^ « z 

positive imaginary axis, the plot may also be used to exhibit the ranges of o 

which, at a given frequency «, correspond to propagating (realgO or non¬ 
propagating (complex df) waves (see (32) et seq.). These propagating wave 
domains have been shaded in Fig. 5; no propagation obtains outside the shaded 
regions. Since the defining expressions for ^ change when ,c^<0, i, e., 
when w<Up (see (A13) ), a separate scale for along the negative vertical 
axis is used in this range. The frequency domain is subdivided naturally into 
the various intervals exhibited in Fig. 5. and denote, respectively, the 
larger and smaller of the frequencies for which s 0, corresponds to 
= 0, while 2£j/(ej + 1) = at 


The behavior of)f^ ^ for waves satisfying the radiation condition at 
z —#+00 in the various frequency ranges is summarized in Table I. The first 
column lists the branch point singularities which lie on the positive real a-axis 
in the ordinary and extraordinary integrals. The notation a^ or a” signifies that 
a branch point lies at a = a, and that the path of integration is indented around it 
into the upper or lower half of the a-plane, respectively. To each a“ there cor¬ 
responds a (-a)^, i. e. , a branch point on the negative real axis, which is not 
listed explicitly. 



PIBMRI-1069-62 


21 


Typical dispersion curves in the various frequency ranges are shown 
in the left half of Fig. 6. These curves represent the real X ■ 0 solutions 
of Eqs. {A9) or (A13) and the portions corresponding to the ordinary and extra¬ 
ordinary waves, as defined in (23a, b), are labeled "o" and "e", respectively.* 

(If no curve is shown, the corresponding wave type does not propagate in this 
frequency range). The singularities of ^^(o) and (o) along the real a-axis 
are exhibited in the right half of the figure. (The branch points at Q- . in (33) 
lie on the real axis only when 6 = - 1) = ^ ^ On the darkened 

segments, X is real; elsewhere it is complex. The proper continuation of the 
functions X (’^) around the various branch point singularities is effected in 

O I 6 

accord with conditions (24) and (31), as summarized in Table I. To render 

these functions unique in the multisheeted complex a-plane, branch cuts have 

been drawn. This partitioning removes any ambiguities as to the disposition 

of the integration path with respect to singularities of X ; the path is dis- 

o, e 

torted around the branch points in a manner illustrated in Fig. 6(a) (a trans- 
formation of variables in (7a, b) from ^ - t| to o can be accomplished via (Bl) ). 
The integrands in (7a, b) will generally also contain pole singularities on the real 
a-axis which must be investigated separately. Such an investigation has been 
carried out in connection with the problem of radiation from an anisotropic half 
space (see reference 8). The discussion herein suffices for the complete analysis 
of radiation in an infinite anisotropic medium where pole singularities do not 
arise. Application of these results to the detailed study of the far field radiated 
by a dipole source is to form the subject of a separate publication (see also 
reference 8). 


The asymptotes of the open branches of the dispersion curves (Figs. 6c, d, f, i) 

can be found by looking for real X, n solutions in (A9) and (A15) as 

Solutions exist only when Cj<0, in which instance \/ ^\ ^ *• 

The open branches arise for the ordinary and extraordinary modes when c >0 and 

z 

4^<0, respectively. 

If the range of propagating modes extends to 0-^± oo (see Figs. 6c, d, f, i), 
absolute convergence of the integrals .can be secured by deforming the end points 
of the integration path away from the real a - axis into a region where the integrands 
decay exponentially. 
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Table I 


Singularities 
(on real cucis) 


Range of propaga¬ 
ting modes 



<o^> «> 


U2> «> Wj 


<i)> max(fc)^<«)^) 


(i) > fa)> <«> 

c p 


rrin(«p, «^) > «> e 


<0 > «> max(w^i w^) 


«•>(<>> u 
4 c 


nin(u^i u^) > u> 0 





o^<l, a^<« 



C^<1, cr^^<c^<oo 


2 2 
a ^<02 


o^<l 


=^<|02l 


a^< 00 
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Appendix A . 

Determination of the vector eigenfunctions 

The transverse electromagnetic fields excited by arbitrary electric 
or magnetic current distributions in the anisotropic medium are to be re^e- 
sented as a superposition of transverse vector modal solutions which indivi¬ 
dually satisfy the source-free field equations and the boundary conditions in 
the transverse domain. To eliminate the z-dependence from (3a, b), charac¬ 
terized by the translation operator a/dz, we assime a separable representa¬ 
tion of the forrn 

- j k if. z » j 

Eti (r) = e. (X. y) e ^ . k = > 0, (Ala) 

- j k3f. z 

(r) = Y. H. (x, y) e (Alb) 


where X. and the normalization parameter are constant, and the subscript 
represents the modal index. Eqs. (3a, b) then define the vector eigenvalue 
problem for the transverse mode functions and 



(AZa) 


AZb) 


Since the transverse (x - y) domain is unbounded, the eigenfunctions have 
a plane wave dependence characterized by 


e. 
—1 




where the continuous spectrum of the normalized wavenumbers 4 and ti runs 
over the values - «> < (^ , n) < oo, and the factor (k/Zir) has been included 


(A3) 
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for convenience in normalization. Since the operator in (A2a, b) is now 
replaceable by -jko_, a = one obtains the algebraic equations 

for the transverse position vectors £ and h (for fixed ^ and q), 

^e = Y4R- (hxz^), YC;«(hx z^) = S • e , (A4) 

where R and S are self-adjoint 2x2 dyadics (note: ~ ®t ^ 

R = 5 ■ it ■ ^-o* ’ (A4a) 

The superscript"^ denotes the adjoint operator^ Eqs. (A4) can be re-expressed 
as the eigenvalue problem 

E • « = e, • (h X z ^) = (h x z ^), P = R • S . 

(A5) 

It follows from the theory of linear operators that the two eigenvectors e and 
il X belonging to the same eigenvalue 7^ satisfy a bi-orthogonality relation as 
in (5). The subscripts o and e are employed to distinguish the two possible eigen¬ 
values, and associated eigenvectors, in (A5). 

Eqs. (A5) are solved readily in a basis comprising the orthogonal 
position vectors £ and £ x , wherein we represent 


£ = a £ + b£ X , h=ca+dcTxz 


(A6) 


A straight-forward calculation leads to 




(1 - ) - j C? (1 - CT > 


c 

d 


= 


c 

d 




(cj - a ) 


Z 12 2 

t a = I + q , 


(A7) 


(A7a) 


from which one obtains the normalized eigenvectors in (4a, b). 
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The eigenvalues 



are the two solutions of the dispersion equation 


) = det (P - V ^) = 0 


(A8) 


or 


^4 . .2 r 2, , , 1 . 4 ,2 2 . . 2 , 2 2 

M H (-1 + 1 ) - 2«1 J + ‘1 - <'1 - '2 + •llo + 'l - '2 

which yields 

TCo = u i - W = u + (ej - i)a/2, 

e 

with U, W, and A defined in {8b) and (4e), respectively. Evidently, 

when U^> W (A real) , 

e e 

and 

^ when < W (A imaginary) . 

e S 



= 0, (A9) 


(AlO) 


(Alla) 


(Allb) 


(A9) is the "Booker quartic" for the longitudinal wave number"^ ( "jf is 

identical with Booker’s q, save for the normalization to k instead of k ; i, e,, 

o 

q= v/e 7l() (see reference 3, Eq. (13.13), with p = 6 = 0). 

The subscripts o and e distinguish the two solutions corresponding 
to the + and - signs in (AlO), respectively. The respective definitions "extra¬ 
ordinary" and "ordinary" are commonly applied in plasma theory to those real 
solutions of the dispersion equation which correspond to waves whose propagation 
in a direction transverse to the z-axis is, or is not, affected by the presence of 
the d. c. magnetic field (along z).* Transverse propagation corresponds to 
= 0; for the ordinary mode one then has a= 1, while for the extraordinary 
mode, a ^ 1. It has been customary to label as "ordinary" and "extraordinary" 
those (real) branches of the'K vs. a (dispersion) curves which do, or do not, pass 
through the points 0» ( 7 = + 1. For our purposes, however, it is more 
significant to effect a definition on the basis of the analytic properties of the 
multivalued functions and (o) as defined in (8a). The resulting appor¬ 

tionment, between and 3^^ , of the real branches of the dispersion curves 
will then not necessarily coincide with that mentioned above (see Sec. IV). 

We have nevertheless retained the terminology "ordinary" and "extraordinary" 

for the o-and e-solutions, respectively, because the real solutions ofdo 

o 

generally include the points = 0 , ± 1. 


* When the magnetic field is absent, e^ = 0 and Ej^ = 1. 
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From (A4) and the orthogonality relations (5), one deduces for 
^he characteristic admittance Y: 


Y 

o, e 



e 



e 


S 


e 

-o» 


e 




o, e 


o, e 


th X7, - R'h xz 
— o, e —o ~ — o, e —o 


(A12) 


Substitution of (4d) into the last expression in (A12) yields formula (8a). 

If one of the regions is isotropic, e, = 1 and £_ = 0 so that £ = 1 . 

2 2 2 ^ ^ 

In this instance. TK- = 3f =l-CJ, and it is convenient to choose the 
'' o e 

conventional linearly polarized E mode (single primes) and H mode (double 
primes) eigenfunctions 

o , . a X z , 

e' = h' X z = - , e" = h" X z = -^ . (A13) 

— — —o o — — —o o 

It can then be shown that the resulting transverse field representation is 

still given by (7a, b) provided that we replace all subscripts o by single primes 

* 

and all subscripts e by double primes. The E and H mode voltages V, V 
and the currents I', I" satisfy the transmission line equations (8), with the 
characteristic impedances defined as 

Z' = = ~~ , Z” = = J1- 0^ . (A14) 

The equivalent network for the mode coupling produced by a plane interface 
between an isotropic and an anisotropic region can then be deduced by proceeding 
as in (10) - (11). 

If £ <0, i. e., k = -j|k|, the normalization of the transverse and 
z 

longitudinal wave numbers q, and is taken with respect to jkj. Hence, 
k in (Ala) and (A3) should be replaced by |k|, whence ^ in (A2a) is defined as 
yj Wq/ I e I • The resulting eigenvalue problem is then found to be the same as 
in (A5) provided that one replaces aby (+ja) and S by (-S). The dispersion 
relation is now given by 

. [cr^(£l + 1) + 2cjj + £jO^ + (£j - £2 + Cj) + £^ - £2 = 0 , (A15) 

i. e., (A9) applies provided that 9(^ and are replaced by (-Jl?) and (-a^),-..^ 
respectively. Finally, the eigenvectors are given as in Eqs. (4), if k and a 

The validity of the formal expressipns obtained via this replacement does 
not necessarily imply that the E and H modes correspond to the limiting 
values of the o- and e- modes, respectively, as £^-*~^l^ . 



PIBMRI-1069-62 A- 

are replaced by Ik( and (+ja), respectively (note: o—^+ja). In summary, 

formulas (7) - (9) appropriate to 0 can be obtained from those for e > 0 

by letting vc a—*jc, 

z z 
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Cylindrical wave representation of the fields 

Instead of the plane wave representation in (7a, b), it is frequently 
more convenient to employ a cylindrical wave representation based on the 
well-known transformation 


I = /de/dr, f(o) * + +jk(ex» +T,y‘) 

.00 -00 


00 . . J 00 

= Zir Yi e * / of (o) J (krrp) J (krrpMdo, k>0, 

n=-oo - n n 


(Bl) 


where 




X = p cos <}), y = p sin <}> , x* = p* cos «|>* , y* = p* sin <])• 

(B2) 


Each of the constituent integrals in {7a, b) resulting upon substitution of {4c, d) 

can be reduced to the form 1 upon observing that the vector o can be replaced 

by the vector operator (j/k) , and that the differentiation and integration 

operations can then be commuted. The scalar constituents of e and h 

— o, e — o, e 

remaimng in the integrands are functions of a only. The voltage and current 
^ and ^ are solutions of the differential equations (8). Since the para¬ 
meters , Y, Z, i, V in this equation depend on {g, q) only through a, as do 
the coupling matrices (11) and (15) descriptive of interface or terminal effects, 
^o, e ^o, e expressed as functions of a only, and a cylindrical wave 

representation is effected as in (Bl). The explicit appearance of the factor 
«XP ^ y ') J fh® integrand of (Bl) results from consideration of a 

point current element located at r* = (x', y', z*) - (p*, <t>«, z*) (cf. (9a, b)). The 

expression for a distributed source is obtained by integration over the source 
point coordinate r'. 

If f(a) = f(-c^, i. e., f is an even function of o, the tr integral can be 
converted to run over the entire real a-axis: 


00 . 

^ f of (d (kap) J^(kap') dor = / of (a) J„(op^) ir^op>) do. 03) 

O 000 


P> ^*^*1 P< denote the greater and lesser of the variables p and p*, respectively. 
The second representation in (B3), involving an infinite integration contour, is 
particiUarly convenient for an asymptotic evaluation of the integral. 
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Appendix C 
Figure Captions 

Fig. 1 - Physical structure and associated network problem 

Fig. 2 - Equivalent network for interface between two gyrotropic media 

Fig. 3 - Definition of traveling and standing wave quantities 

Fig. 4 - Dispersion curve, wave normal, and ray 

Fig. 5 - Frequency dependence of various plasma parameters 

Fig. 6 - Dispersion curves and singularities on Re o-axis 

The dispersion curves on the left show the behavior of vs. O 
for propagating waves in various frequency ranges; the designations 
"o" and "e" distinguish the ordinary and extraordinary branches, 
respectively. Singularities due to 9 ^ (o) on the real ft- axis are 

shown on the right. The dark segments denote propagating wave 
regions. 


THE PHYSIGAL STRUCTURE AND THE ASSOCIATED NETWORK PROBLEM 
(LONGITUOINAL GYROTROPC AXIS ) 
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